Anisotropic states of two-dimensional electrons in high magnetic fields 
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We study the collective states formed by two-dimensional electrons in Landau levels of index 
n > 2 near half-filling. By numerically solving the self-consistent Hartree-Fock (HF) equations for a 
set of oblique two-dimensional lattices, we find that the stripe state is an anisotropic Wigner crystal 
(AWC), and determine its precise structure for varying values of the filling factor. Calculating the 
elastic energy, we find that the shear modulus of the AWC is small but finite (nonzero) within the 
HF approximation. This implies, in particular, that the long-wavelength magnetophonon mode in 
the stripe state vanishes like q 3 ^ 2 as in an ordinary Wigner crystal, and not like q 5 ^ 2 as was found 
in previous studies where the energy of shear deformations was neglected. 

PACS numbers: 73.43.Nq, 73.20.Qt 



One of the most fascinating properties of two- 
dimensional electrons in a strong perpendicular mag- 
netic field is the strong anisotropy these systems display 
around half filling of Landau levels (LL's) of index n > 2. 
This anisotropy, first seen in DC transport measurements 
0, Q, has been interpreted early on 0, 0] as evidence 
for the existence of a unidirectional charge density wave, 
termed "stripe state" , near half filling, whereas the insu- 
lating states observed for small current at intermediate 
and small partial filling factor v* = v — 2n iv being the 
total filling factor) were identified as pinned crystalline 
states of the 2D electron system HUES 00- While 
the concept of a stripe phase, with the density of elec- 
tron guiding centers taking alternating uniform values 
that are piecewise constant along a periodic stripe struc- 
ture, is sufficient to explain the anisotropy of the DC 
response, other structures may lead to an anisotropic 
response as well. One such structure is an anisotropic 
Wigner crystal (AWC), whose existence was suggested 
by microscopic Q as well as Hartree-Fock (HF) 0, E3 
studies, and which under certain circumstances (see be- 
low) resembles a stripe structure, but with a periodic 
modulation of the electronic density along the stripes. 
In this Letter, we report on a systematic study of the co- 
hesive energy of anisotropic triangular WCs, with special 
focus on the LL n — 2. By numerically solving the self- 
consistent HF equations, we find that the stripe state 
is in fact an AWC. Indeed, as it can be seen in Fig. 



^ a highly anisotropic rhombic lattice may be viewed 
as a periodic arrangement of one-dimensional stacks of 
equidistant electron guiding centers (see Fig. P). Calcu- 
lating the elastic moduli of the AWC, we find that the 
energy of shear deformations corresponding to channels 
sliding past each other, although very small, does not 
vanish. This has the important consequence that the 
stripe state behaves as a smectic on short length scales, 
but on scales longer than a crossover length L cr (that we 
estimate below, see Eq.©) it behaves like an ordinary 
two dimensional anisotropic solid. 

In order to probe anisotropic crystalline states and 
the corresponding cohesive energies of the 2D electron 
system, we have used the self-consistent Hartree-Fock 
method of Ref. In this work, we have system- 

atically studied a family of oblique (rhombic) crystals 
parametrized by the following lattice vectors (n and m 
are integers): 

R™ = nai x + (2m + n)a 2 y , (1) 

where a\ — \[Zaj2\[\^~e and a 2 = \J\ — ea/2, and 
where e is a positive parameter (such that < e < 1) 
which quantifies the degree of anisotropy of the lattice 
(with e = corresponding to an isotropic triangular lat- 
tice), and a = £(AttM/V3 v*) 1/2 is the average lattice 
spacing of an isotropic lattice with e = at the same den- 
sity (M is the number of electrons per bubble in a general 



2 



2 ....... 2 



....... . o 



_2 ....... . _2 



-6 -4 


-2 

e = 


2 4 6 -6 -4 -2 




e = 


2 4 

0.25 


4 

2 

2 
-4 






















2 
O 
-2 

















-6 -4 -2 O 2 4 6 -6 -4 -2 2 4 6 



e = 0.5 e = 0.75 

FIG. 1: (Color online) Evolution of the WC appearance for 
various values of the anisotropy parameter e of Eq. Q. For 
e > 0.5, the lattice shows a pronounced channel-like structure, 
with each channel consisting of a one-dimensional periodic 
chain of electron guiding centers. 

bubble crystal). At any fixed partial filling factor v* , we 
compute the cohesive energy E co h{e) for a sampling of 
values of e in the interval < e < 1, and determine the 
value of £ that minimizes this last quantity. We find that 
the cohesive energy of bubble crystals with M = 2 and 
3 is minimized by the value e = throughout the whole 
v* range. For the Wigner crystal (with M = 1) however, 
we find that the cohesive energy is minimized by nonzero 
values of e for v* > 0.22. The resulting variation e(v*) 
for the WC is shown in the upper panel of Fig. |3 while 
the lower panel shows the resulting cohesive energies of 
the various phases. As it can be seen, the ground state 
of the 2D electron system goes through a succession of 
isotropic ground states (Wigner crystal and bubble solids 
with M = 2 and 3) at small and intermediate values of 
the partial filling factor v* . At v* ~ 0.42, the system 
undergoes a first order phase transition to a modulated 
stripe state, which in this case is nothing but an AWC, 
with values of e ranging from e ~ 0.76 at the transition 
to about e ~ 0.80 at half filling in LL n = 2. 

Having investigated the phase diagram of the 2D elec- 
tronic system, we now want to investigate the elastic 
properties of the resulting crystalline phases. For the 
anisotropic triangular lattice we are considering here, the 
general theory of linear elasticity in two dimensions dic- 
tates that the elastic energy density be of the form |Tl| : 

S e i = i (cnu^ + c 22 U22 + 4c 66 ui2 + 2C12U11U22) , (2) 

where u a p — ^(d a up + df3U a ) is the symmetric strain ten- 
sor, and where we have used reflection symmetry (in both 
x and y) to reduce the number of elastic moduli from six 
to four. (For an isotropic triangular crystal, the number 
of elastic moduli can further be reduced by symmetry 
[llj . so that the elastic energy can be written in terms 
of only two elastic constants, the compression (en) and 
shear (cqq) moduli.) Note that stability of two dimen- 
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FIG. 2: (Color online) Upper panel: values of the anisotropy 
parameter e that minimize the cohesive energy of the WC 
vs. v* . Lower panel: cohesive energies of the Wigner crystal 
(circles), 2e bubble solid (squares), and 3e bubble solid (tri- 
angles). For the WC, e values from the upper panel are used. 
Both panels are for LL n = 2. 



sional solids requires that the shear modulus be strictly 
positive, C66 > 0. For two-dimensional uniform deforma- 
tions, such that the displacement vector u(r) is given by 
u a = Un.nXg, with constant coefficients u a ^, it can be 
shown |l2( that the reciprocal lattice also experiences a 
homogeneous deformation, with the deformed reciprocal 
lattice vectors {Q} given in terms of the original ones 
{Q} by Q a = Q a — up.aQp + 0(u 2 ). To extract a given 
elastic constant in Eq. we calculate the HF energy 
of the corresponding distorted state (with distortion am- 
plitude ito ), and extract the elastic constant from the 
excess (elastic) energy AE e i(uo) = Ehf(uo) — Ehf(0). 
For example, to extract the shear modulus cqq we cal- 
culate the HF energy Ehf(uq) of the distorted crystal 
(which for small uq has an optimal anisotropy e that is 
very close to that of the undistorted lattice) using the 
reciprocal lattice vectors (Q x = Q x — uqQ v ,Q y = Q y ), 
corresponding to the shear deformation polarized along 
y such that u(r) = UQxy. The shear modulus (in units 
of energy /particle) is then given by cq§ — 2AE e i(uo)/ug 
(to obtain the shear moduli in the conventional units of 
energy per unit area, one should divide by the area of 
a unit lattice cell A c = 2ttM£ 2 /v*). In our case where 
we consider electrons interacting through the long range 
Coulomb potential, the above procedure is appropriate 
only for the shear moduli (which can be shown to be only 
weakly dispersive 0, and for the non-dispersive 

part of the compression moduli (the compression mod- 
uli are well-known to be dominated by a l/<; plasmon 
contribution in the long wavelength limit). 

We now turn our attention to the calculation of the 



3 



shear modulus of the actual anisotropic ground states of 
our system. We choose to do the calculation in two differ- 
ent ways, by calculating the energy of shear deformations 
both along x (u(r) = u yx) and along y (u(r) = u$xy). 
In Fig. |31 we plot the behavior of the two correspond- 
ing moduli, that we call cqg^ x and cee tV vs. v* for the 
AWC ground states we found above in the n = 2 LL. We 
first note that cgq iX and cee iV are equal, as they should 
be, which confirms the fact that the shear deformations 
in our lattice are controlled by a single shear modulus 
C66- Comparing the values of this shear modulus Cqq ob- 
tained over the stripe phase to the values obtained in 
Refs. ElQfor the Wigner crystal and bubble phases, we 
see that cqq jumps discontinuously by over two orders of 
magnitude from values of order 0.01 (e 2 / nt) for v* < 0.42 
where the electronic solid is isotropic, to values of order 
5 x I0~ 5 (e 2 / nl) in the anisotropic regime 0.42 < v* < 0.5 
(k is the dielectric constant of the host semiconductor). 
This indicates that sliding the channels in Fig. past 
each other costs a finite, albeit very small shear energy. 
One may of course worry that the smallness of the shear 
moduli we find above is a numerical artifact with no ac- 
tual physical meaning. In the following, we shall argue 
that non-zero shear moduli are a natural consequence of 
the modulated density along the stripes in the AWC. 

Central to our discussion is the observation that there 
is a fundamental difference between our AWC state, 
where electron guiding centers form modulated channels 
breaking translational invariance in both x and y direc- 
tions, and the conventional stripe state, with a uniform 
density of guiding centers along the stripes. For, while 
in the latter case, assuming the stripes are parallel to 
the y direction, a deformation of the form u = Uoyx 
corresponds to a rigid rotation of the whole stripe struc- 
ture (with no overall energy cost), in our case with the 
electron guiding centers forming modulated channels, the 
deformation u = itoyx does not correspond to a pure 
rotation of the system, and in fact involves a compres- 
sion of the electron guiding centers along the direction 
of the channels, which costs a nonzero elastic energy 
~ cee(d y u x ) 2 . This point can be seen more clearly if 
we consider the following displacement field of the AWC, 
Unm = —0(n + 2m)a2X which, in continuous notation, 
corresponds to the shear deformation u(r) = —9yx. 
Performing an infinitesimal rotation of angle 9 to the 
new cartesian frame (x'y') such that x = x' — By' and 
y = Ox.' + y', we obtain that the displaced positions 
r nm = R nm + u„ m are given in the {x'y') frame by 
Trim = naix' + [(n + 2m)a,2 — na\0\ y' . As it can be seen, 
the shear deformation u nm = — 9(n + 2m)a2X introduces 
a compression of the electron guiding centers along the 
direction of the channels (i.e. along y'), hence the finite 
energy cost of such a deformation. By contrast, a rigid 
rotation with displacement vectors of the form u nm = 
9\na\y — (n + 2m)a2x] leads to displaced lattice vectors 
(to order 9 2 ) of the form: r nm — nctix' + (n + 2m)a,2y' ', 
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FIG. 3: (Color online) Shear modulus of the anisotropic 
Wigner crystal in n = 2 in units of e 2 / nl. The values of 
C66 obtained for a shear polarized along x and along y are 
denoted by cee,x and cm,y, respectively. 



which corresponds to a perfectly ordered anisotropic WC 
in the (x'y') frame. The above analysis clearly shows 
that, while for a conventional liquid crystal (or stripe) 
with a uniform density along the stripes a deformation 
of the form u(r) = —Oyx corresponds to a rigid rotation 
of the system, in our case where translational invariance 
is also broken along the direction of the channels the 
above displacement field does not correspond to a rigid 
rotation, and in fact introduces a compression of the elec- 
tron guiding centers along the stripes with a finite elastic 
energy cost. 

The elastic matrix of the AWC has elements given by 
^(q) = cu{q)q 2 x + c m q 2 y + Kq*, $ OT (q) = c 22 {q)q 2 + 
c 6 6<?z, and <& xy (q) = ®yx(q) = (012(g) + c^)q x q y , where 
we have introduced a term ~ ^Qy i n to account for 
the bending rigidity of the stripes at small length scales. 
In conventional smectics, where the shear modulus cqq is 
exactly zero, the term ~ Cee(d y u x ) 2 in the elastic free en- 
ergy is replaced by the bending rigidity term ~ K(d 2 u x ) 2 
[l5j. In our case, where Cqq is very small but nonzero, the 



elastic energy of the AWC will contain both terms, and 
phenomena taking place on various length scales will be 
governed by one or the other of these two elastic con- 
tributions. In particular, one can define a characteristic 
wavevector q cr = y^cee/K (or, equivalently, a crossover 
length scale L cr ~ l/q cr ) which delimits the long wave- 
lengths regime (0 < q -C q c ) where c^q 2 ^> Kq 4 and 
where the elastic properties of the system are governed 
by the ordinary elastic energy of a two-dimensional lat- 
tice, and the short wavelengths regime q c < q < qsz 
where the elastic properties of the system are smectic- 
like. Using the value K = 189 mK obtained in Ref. 
and taking cqq ~ 5 x 10~ 5 (e 2 /e£)/A c , we obtain the fol- 
lowing estimate for the crossover length scale: 

L cr ~10a. (3) 

As an illustration of the qualitative argument above, 
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in the following we shall discuss the behavior of the nor- 
mal modes taking into account the shear modulus cqq of 
the AWC. In the presence of a magnetic field, the equa- 
tion of motion for the electron at lattice site can be 
written in the form (xij is the displacement vector of the 
electron with respect to its equilibrium position Rj, m* 
is the effective mass of the electron in the host semicon- 
ductor and e a p is the totally antisymmetric tensor in two 
dimensions) 



m 



, d 2 u tc 



j 



eB du il3 



We shall seek a solution to the above equation of mo- 
tion that represents a wave with angular frequency u> and 
wavevector q, i.e. u a (Rj) = A a (q)e^ q -' Ri ^ UJt \ where the 
A Q 's are complex coefficients whose ratios specify the rel- 
ative amplitude and phase of the vibrations of the elec- 
trons within each primitive cell. Substituting the above 
expression into Eq. (0J results in the following secular 
equation for the Wigner crystal: 



D xx (q) 



D xy (q) 



D xy (q) + iLULU c D yy (q) 



10JLU C 
, ,2 



0, (5) 



where D a p is the dynamical matrix, which is defined 
in terms of the elastic matrix <I> Q) g by D Q ^(q) = 
-^y\$ Q/3 (R,) e -^.. 

The above secular equation has a non-vanishing solu- 
tion (for the A Q 's) only if the determinant of the secular 
matrix is zero. This leads to the following eigenmodes: 



: (q) 



D xx + D yy ± VA(q) 



(6) 



-D„ V ) 2 +4D 2 



where A(q) = ul+2ul(D xx +D y y)+(D xx ~^ ml , 
The solutions with a (— ) sign represent magnetophonon 
modes, with eigenfrequencies which vanish as q — > 0, 
while the solutions with a (+) sign correspond to mag- 
nctoplasmon modes which tend to oj c as q — ► 0. Solving 
explicitly for the magnetophonon mode, we obtain in the 
long wavelength limit where c\±(q) ~ 022(9) — c i2(<?) = 
e(q) the result ul (q) = ^ [q x - 2q 2 x q 2 y + q §. Given 
that cn(q) ~ l/g, this implies that magnetophonon 
eigenmodes have the characteristic Wigner crystal behav- 
ior to(qi) — w(qy) ~ q 3 ^ 2 for q < q cr . On shorter length 
scales, the bending rigidity term dominates over the shear 
term. Neglecting c 66 in the expression of the elastic ma- 
trix, we obtain w?_(q) ~ {Kcu(q) /m 2 uj 2 )qy, with the 
consequence that m-(qy) ~ q 5 ^ 2 for q > q cr . In order 
to test these predictions, we have performed a direct nu- 
merical calculation of the dispersion curve using the gen- 
eralized random phase approximation (GRPA) method 
ofRef. H The results of this computation are shown 
in Fig. 0] and seem to indeed confirm the results of the 
above elementary calculation (although the best fit we 
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FIG. 4: (Color online) Dispersion curve for the stripe state 
at filling factor v = 4.45. There is a change of behavior from 




q 1 - 8 to lj 



q ' around q cr = l/10a (in the y axis label, 



lj is expressed in unit such that h = 1). 



obtain at small q for a law of the form to ~ q a is with an 
exponent a ~ 1.8 instead of 1.5). 

In conclusion, in this paper we have investigated the 
possibility that the stripe state near half filling of higher 
Landau levels (with index n > 2) may be anisotropic 
crystalline states. By numerically solving the HF equa- 
tions for a family of oblique realizations of the electronic 
lattice parameterized by an anisotropy parameter e, with 
< e < 1, and finding the values of e that minimize 
the cohesive energy of the system, we find within a self- 
consistent HF approach that the ground state of the sys- 
tem in the stripe phase is indeed an anisotropic Wigner 
crystal, with a small but finite shear modulus. As an 
immediate consequence of this, we have shown that the 
magnetophonon modes have the usual dispersion of an or- 
dinary Wigner crystal w(q) ~ q 3 ^ 2 at wavelengths longer 
than a crossover length scale L cr , which we estimate to 
be of order 10 lattice spacings near half filling in n = 2, 
with w(q) ~ q^l 2 on shorter length scales. Previous scale- 
sensitive predictions for the stripe state, such as predic- 
tions of the electromagnetic response of the stripe state 
in presence of pinning disorder, will only hold at length 
scales shorter than L cri with the standard Wigner crystal 
behavior taking over at longer length scales. 
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